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Calculation  of  Fission  Neutron  Age  in  NaZrF^ 


One  of  the  limits  imposed  on  experimental  measurements  of  the  neutron 
age  of  various  materials  is  the  size  required  for  the  experimental  setup  if 
the  age  is  too  great.  In  the  case  of  NaZrF5,  a  calculation  has  been  made 
to  determine  if  this  is  a  limiting  factor  in  the  measurement  of  age  to 
indium  resonance.*  In  the  calculation  it  is  assumed  that  it  is  possible  to 
apply  the  Fermi  age  theory  for  the  distribution  of  neutrons. 

Method  of  Calculation 

For  a  monoenergetic  point  source  of  neutrons^  the  solution  of  the  age 


equation  is^ 


where 


,(r,Tr)  .  a 


r  =  distance  from  point  source, 

Q  =  nuinber  of  neutrons  emitted  per  second  by  the  source, 
q  =  neutron  slowing  down  density. 

'The  age  ^  is  defined^  by 


'C(Eo,E)  = 


3  ^ (l  -  COS&) 


where 


E  =  neutron  energy, 

Eq  =  source  energy, 

X(e)  =  mean  free  path  at  energy  E, 

^  =  average  loss  in  logarithm  of  energy  upon  collision, 

COS0  =  average  cosine  of  angle  between  neutron  direction  before  and  after 
collision  in  laboratory  system. 


*  The  age  to  indixam  resonance  is  the  square  of  what  is  usually  called  the 
"slowing  down  length”  to  indium  resonance. 

1  Jay  Qrear,  A.  H.  Rosenfeld,  and  R.  A.  Schluter,  Nuclear  Physics,  IMiversity 
of  Chicago  Press,  Chicago,  1950- 
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The  mean  square  slovdjng  down  length,  p,  is  given  by 

r 

—r  q(r,r)dV 


=  10 


(3) 


't)dV 


/'  Q  r, 

/  I^«?)3/2  ®  4«r2dr 


/ 


'o 


(3a) 


dr 


In  the  case  In  which  the  source  Is  not  monoenergetic,  its  energy  distribution 
being  described  by  f(E), 


q(r,E) 


j  dE'f(E') 
E 


Q  e-rV  4|t(E’,E")-'t(E,E'')] 
(4k  j2(E',E")  -  T(E,E")~|1 


(4) 


where  E"  is  some  lower  energy.  If  r2(e)  ie  computed  from  Eqs.  (3)  and  (4) 
it  is  found  that 


R^(e) 

6 


OD 

J  fCEOdE'tCE^E) 

E 


(5) 


The  right-hand  expression  of  Eq.  (5)  will  be  called  the  age  to  energy  E  for 
the  distribution  f(E).  Integrating  by  parts. 


i 


f(E')dE'  t(E',E)  =  -  t(E',E)  /  f(E")dE" 


dE'  A-  't(E',E)  /  f(E")aE’’ 

'  dE'  J 

E  E' 


It  is  reasonable  to  assume 
Eq.  (6)  vanishes.  Thus 


that  the  first  term  on  the  right-hand  side  of 


't(E)  =  /  dE’  iX  (E’,E)  /  f{E’*)dE”  W) 

/  dE’  7 

E  E’ 

where  "^Ce)  is  the  age  to  energy  E  for  the  distribution  f(E).  If  the  value 
of  dV*^  is  derived  from  Eq,.  (2)  and  substituted  in  Eq^.  (?)> 

/oo 

A')dB’  /  f(E'/)dE"  =  ^(E)  (8) 

5fE’(l  -  COS0)  J 
E  '  E' 

For  the  case  of  a  fission  spectrum 


f(E)  =  -x /—  e**^  sinh  'V2E 

V  Tte 


where  E  is  the  neutron  energy  (Mev).  It  is  easier  for  the  purposes  of 
integration  to  change  variables  by  the  transformation  E  =  V^/2.  Thus 


f(E)dE  =  e  sinhV  VdV 

tJte 


-5- 


-6- 


The  advantage  of  Eq..  (ll)-)  is  that  G(V)  may  be  expressed  analytically 
while  the  other  terms  may  be  directly  computed  from  the  cross  sections.  This 
with  Eq,  (14)  it  would  only  be  necessary  to  do  one  numerical  integration.  With 
Eq,  (5)  it  would  be  necessary  to  first  compute  ^(E^^E)  by  numerical  integration 
and  then  perform  a  second  numerical  integration  to  obtain  X(s)*  ^ 

other  hand^  it  is  sometimes  true  that  ^  (E’,E)  is  available  from  experi¬ 
mental  data,  in  which  case  it  would  be  advantageous  to  use  Eq,  (5)* 

Numerical  Results 

Values  of  G(v)  that  were  computed  are  given  in  1?able  1,  To  aid  in  the 
integration,  polynomial  fits  were  made  to  G(v). 


G(V) 

=  2.5066 

+  0.1592V  -  0.4708v^ 

for  0 

V  ^1, 

(15) 

0 

E  —0.5 

G(V) 

=  3.3535 

-  1.1785V 

for  1  a 

V  i.2 

(16) 

0.5  £ 

E  £.2 

G(V) 

=  2.6053 

-  0.8044V 

for  2  ^ 

V  £3 

(17) 

2  :S 

E  i.4.5 

G(V) 

=  0.7249 

-  0.1776V 

for  5  ^ 

V  S.  4 

(18) 

4.5  £ 

E  8 

G(V) 

=  0.0709 

-  0.0141V 

for  4  i 

Vi:  5 

8  - 

E;l2.5 

(19) 

Equations  (15)  through  (19)  give  polynomial  fits  to  G(v)  for  the  ranges  of  V 
indicated.  The  corresponding  values  of  the  energy  in  Mev  are  also  indicated. 
The  error  in  the  quadratic  fit  is  about  Ijt.  The  maximum  error  in  Eq.  (I6)  is 
about  lO^t.  In  the  remainder  of  the  linear  fits,  the  maximum  percentage  errors 
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Table  1 

Values  of  G(v) 


V 

G(V) 

V 

G(V) 

0.0 

2.5066 

1.0 

2.1750 

0.1 

2.5062 

1.5 

1.4328 

0.2 

2.5034 

2.0 

0.9965 

0.3 

2.4959 

2.5 

0.4905 

0.4 

2.4815 

3.0 

0.1921 

0.5 

2.4585 

3.5 

0.0595 

0.6 

2.4253 

4.0 

0.0145 

0.7 

2.3808 

4.5 

0.0028 

0.8 

2.3243 

5.0 

0.0004 

0.9 

2.2556 
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become  steadily  worse,  but  the  absolute  error  is  about  the  same.  Further¬ 
more,  in  the  integration, G(V)  is  weighted  by  a  factor  of  l/V  so  that  the 
lower  values  of  V  are  more  important.  The  value  of  Eq  that  will  be  used  is 
1.44  ev,  the  indium  resonance  energy.  The  corresponding  value  of  V  is  I.697  x 
10"^.  Strictly,  the  normalizing  factor  in  Eq,.  (l4)  should  be  G 
instead  of  l/ST,  but  since  Eq  is  so  low  no  appreciable  error  is  introduced. 

To  evaluate  Eq.  (l4)  with  the  polynomial  fits,  another  approximation  was 


made.  By  the  theorem  of  mean  value 


dV(2X^/v)G(V)  ^  _ 

3  ^(1  -  cos?)  ySiT 


- _ 

_ _  f  GMdv  (20) 

1  -  cos0)i&jr  j  V 

J  Vi 


where  the  har  over  the  expression  in  "brackets  denotes  some  value  between  and 
V2»  If  the  polynomial  fits  are  used,  the  integral  on  the  right-hand  side 
of  Eq.  (20)  may  be  evaluated  exactly.  It  is  then  possible  to  obtain  upper 
and  lower  limits  for  the  integral  on  the  left-hand  side  by  using  upper  and 
lower  limits  for  the  expression  in  brackets. 

To  test  the  accuracy  of  this  method,  the  age  in  graphite  (density  = 

1.6  g/cc)  was  calculated.  The  quantities  f  and  cos?  are  constant.  Table  2 
gives  the  maximum  and  minimum  values  of  the  quantity  I  from  which 


the  age  for  graphite  can  be  determined: 


(1.31) 


/  I  9N0  X  10"  J  I  3  5  (1  -  COS0) 


i0.9k) 


J  I  pNo  X  lO’^M  1  3  ^(1  -  COS0) 


=  362  cm^  (max,) 


=  260  cm^  (min.) 


where 


A  =  atomic  weight  of  carbon. 
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Table  2 


Maximum  and 


Minimum  Values  of 
for  Graphite 


E(Mev) 

f  G{Y)dV 

X,2  X  102 

(barn“^) 

Max. 

Min. 

J  V 

X  10"^  -  .08 

15.710 

4.75 

4.73 

o 

• 

o 

00 

1 

o 

Cn 

1-947 

8.5 

.  4.73 

0.5  -  2 

1.146 

34.6 

8.5 

2  -  oo 

0.285 

34.6 

34.6 

SuTii  nf  /?(V)dV  ^  ^2 

1.31 

0.94 

^  V 

p  =  density  of  carbon  =  1.6  g/cc, 

Nq  =  AVQgadro*s  number. 

o 

The  average  between  the  upper  and  lower  limits  gives  an  age  of  311  cm  for 

3  2 

carbon  which  differs  by  only  6^  from  the  experimental  value  of  330  cm  . 
Application  of  Method  to  NaZrFc 

In  an  application  of  this  method  to  the  compound  NaZrF^  the  density  was 
taken  as  h  g/cc/  and  the  quantities  and  cos0  were  taken  as  the  average  of 
the  quantities  for  Na,  Zir,  and  F  weighted  hy  the  cross  sections  and  the  relative 
amount  of  each  element  in  the  compound.  Table  5  gives  the  maximum  value  of 
the  quantity  x^/j  ^(l  -  cos?)  for  the  energy  range  listed;  Table  1;  gives  the 
minimum  value  for  the  same  quantity.  In  these  calculations  12.5  ='0°.  From 


Table  5  the  upper  estimate  on  the  age  may  be  obtained  by 


(89.572  X  10’^) 


where 


A  =  molecular  weight  of  NaZrF^, 

^  s=  density  of  NaZrF^  =  k  g/cc, 

Nq  =  Avogjadro’e  number. 

Similarly,  the  lower  limit  is  estimated  from  Table  4  as  328  cm^.  The  average 
between  the  upper  and  lower  limits  is  435  cm^  which  gives  an  estimate  of  the 


age  in  NaZrF5  as 

r(l.i|^  ev)  =  h33  +  105  =  435  (l  +  O.25)  cm^  (2l) 

The  number  105  is  the  difference  between  the  mean  and  either  the  upper  or 
lower  estimates. 


5  J.  E.  Hill,  L.  D.  Roberts,  and  G.  E.  McCammon,  "Slowing  Down  of  Fission 
Neutrons  in  Graphite,"  ORNL-I87  (Jan.  19^  19^9)* 

4  W.  K.  Ergen,  personal  communication. 
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Table  5 

Maximum  Value  of 

J  3/(1  -  5os0)V 

for  NaZrF5 


E(Mev) 

^^/5f(l  -  COS0)  X  105 
(betrn”^) 

% 

t  'J(V) 

{  V 

1 

1.44  X  10""^  to  2  X  10~5 

5.5 

9.0770 

2  X  10"5  to  2.5  X  10"^ 

2.8 

0.2803 

2.5  X  10-5  to  5  X  10"^ 

0.9 

0.2293 

3  X  10“^  to  4  X  10"5 

2.5 

0.3616 

4  X  10“^  to  6  X  10’^ 

3.0 

0.5099 

6  X  10-5  to  2.5  X  10-2 

4.4 

1.7957 

2.5  X  10"^  to  1.0  X  10’^ 

5.3 

1.7334 

1.0  X  10-^  to  1.1  X  10"^ 

0.5 

0.1176 

1.1  X  10“^  to  2.5  X  10“^ 

3.5 

0.9962 

2.5  X  10"^  to  5.5  X  10’^ 

1.7 

0.3928 

3.5  X  10"^  to  5  X  10'^ 

2.4 

0.5390 

5  X  10"^  to  2.0 

13.0 

1.1460 

Minimum  Value  of 


Table  4 

f  G(v)\2dV 


3f(l  -  cose)V 


for  NaZrFi 


x2/3f(l  -  cose)  X  10- 
(barn"^ ) 


2.0 

‘^.5 


■  X  10"^ 

to 

1,6  X  10"^ 

X  10"^ 

to 

2.2  X  10"5 

X  10“^ 

to 

4  X  10'5 

10“^ 

to 

6  X  10‘5 

10-5 

to 

1.5  X  10'^ 

X  10"^ 

to 

2.4  X  10"^ 

X  10"^ 

to 

8.3  X  10"' 

X  10"^ 

to 

1.35  X  10"^ 

>  X  10”^ 

to 

2.1  X  10'^ 

X  10'^ 

to 

4  X  10"^ 

10“1 

to 

5  X  10"^ 

10’^ 

to 

2.0 

to 

4.5 

to 

oo 

Sum  of  -  cose)  x  10^ 


8.7965 

0.^^001 

0.71^87 

0.5099 

1.1551 

0.5912 

1.5555 

0.6007 

0.14-791 

0.7525 

0.24714 

1.1460 

0.2519 

0.0327 


all),  avi 

V 


54.522 
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